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We investigate a phenomenological model for the spin glass phase of La2- x Sr x Cu04, in which it 
is assumed that holes doped into the CuC>2 planes localize near their Sr dopant, where they cause a 
dipolar frustration of the antiferromagnetic environment. In absence of long-range antiferromagnetic 
order, the spin system can reduce frustration, and also its free energy, by forming a state with an 
ordered orientation of the dipole moments, which leads to the appearance of spiral spin correlations. 
To investigate this model, a non-linear sigma model is used in which disorder is introduced via 
a randomly fluctuating gauge field. A renormalization group study shows that the collinear fixed 
point of the model is destroyed through the disorder and that the disorder coupling leads to an 
additive renormalization of the order parameter stiffness. Further, the stability of the spiral state 
against the formation of topological defects is investigated with the use of the replica trick. A critical 
disorder strength is found beyond which topological defects proliferate. Comparing our results with 
experimental data, it is found that for a hole density x > 0.02, i.e. in the entire spin glass regime, the 
disorder strength exceeds the critical threshold. In addition, some experiments are proposed in order 
to distinguish if the incommensurabilities observed in neutron scattering experiments correspond to 
a diagonal stripe or a spiral phase. 

PACS numbers: 75.10.Nr,74.72.Dn,75.50.Ee 



I. INTRODUCTION 



Generalities 



This paper discusses the influence of disorder on the 
properties of weakly hole doped cuprate materials. In 
cuprates, the superconducting state emerges through 
chemical doping of a parent compound which is insu- 
lating and shows antiferromagnetic (AF) order with a 
high critical Neel temperature of typically a few hundred 
Kelvin. As a consequence of chemical doping, the com- 
pounds are intrinsically disordered. Especially at weak 
doping concentrations, disorder is known to strongly in- 
fluence the behavior of these materials. This is evident 
in the simplest cuprate superconductor, La2_ x Sr x Cu04, 
where the superconducting phase emerges via doping di- 
rectly from a low temperature spin glass (SG) phase. Re- 
cently, glassy characteristics were detected even inside 
the superconducting phase (see Ref. 1 for a summary of 
the available experimental data). 

Understanding the very weak doping regime of 
cuprates, the insulating AF and SG regime, should be 
relatively simple. This optimism is based on the belief 
that this regime is dominated by the behavior of iso- 
lated holes in presence of well developed AF moments. 
The single hole properties seem now to be quite well un- 
derstood and early theories of high temperature super- 
conductivity were constructed from these one-hole wave 
functions. Shraiman and Siggia^ proposed a theory of 
interacting hole-quasiparticles based on the one-hole pic- 
ture and predicted the formation of spiral correlations 
with a pitch proportional to the hole density. Exper- 



iments have to date however not found any evidence of 
such spiral correlations inside the superconducting phase. 
The pairing mechanism suggested by this semiclassical 
picture, a dipole-dipole interaction between holes medi- 
ated by soft spin waves^ 4 ^ has, perhaps unfairly, received 
scant attention of late. A potential weakness of the ap- 
proach is the semiclassical treatment of spins (large S), 
which implies the assumption of a large AF correlation 
length, whereas in the superconducting phase the spins 
are believed to form some kind of quantum disordered 
spin liquid. The scattering of holes by spin excitations 
would then be qualitatively different at large scales. How- 
ever, while the semiclassical theory is formulated for large 
scales, the structure and energy of the resulting two- 
hole bound state is determined by the shortest cutoff in 
the system^, where AF correlations are still intact. Fur- 
thermore, the correlation length can be substantial even 
in superconducting samples, e.g. it exceeds 200 A in 
the stripe compound Lai.45Ndo.4Sro.i5Cu04i» Thus, the 
pairing mechanism suggested by the semiclassical picture 
may hide some truth despite the absence of long range 
order. 

While a semiclassical approach to the superconduct- 
ing regime may or may not be valid, at sufficiently low 
hole concentrations, where static AF correlations are still 
dominant, i. e. in the SG and AF phase, a semiclassi- 
cal treatment of spins is certainly justified. However, at 
these low densities, where the system is still a Mott insu- 
lator, screening is very poor and long-range Coulomb in- 
teraction leads to a strong disorder potential which must 
be taken into account. Here we discuss a model in which 
the entire charge distribution is assumed to be quenched. 
Each hole, localized close to an ionized dopant, is as- 
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sumed to produce a long-ranged dipolar-shaped frustra- 
tion of the AF, similar to the one known to be produced 
by delocalized holes. A polarization of the dipole mo- 
ments then implies the appearance of spiral correlations. 

It is known that the spiral state described by Shraiman 
and Siggia, if one ignores disorder, is unstable toward 
a local enhancement of the spiral pitch. This instabil- 
ity arises from the fermionic susceptibility of the holes 
and may signal an instability towards charge density for- 
mation or phase separation^ However, if the holes are 
quenched this instability is suppressed. Therefore, dis- 
order takes a prominent role in the creation of a spiral 
state. 

We here develop a renormalization approach for dis- 
ordered spiral phases, where we study the scaling of the 
spin stiffness and of the disorder. The importance of 
topological defects of the spiral texture is analyzed and 
their relevance for the physics of the spin glass phase is 
discussed. 



B. Undoped and weakly doped cuprates 

Undoped La 2 Cu0 4 is a charge transfer insulator with 
an antiferromagnetically ordered ground-state. It is well 
described by a simple square lattice spin- 1/2 Heisenberg 
model, 

H B = J^8i'8 jt (1) 

<«) 

with the antiferromagnetic exchange J ~ 1200K. The 
sum is over nearest neighbor pairs of sites and Si are 
spin- 1/2 operators. 

In the study of magnetism of La2Cu04, an approach 
based on the quantum-non-linear-er model (QNLoM) has 
been highly successful. It correctly describes the long 
wavelength hydrodynamic modes (spin waves) of the 
Heisenberg models In this continuum model, it is as- 
sumed that the antiferromagnetic correlation length is 
much larger than the lattice spacing and the model de- 
scribes slow fluctuations of the locally well defined stag- 
gered magnetization n (with n 2 = 1). The QNLerM ac- 
tion is 

¥=tf ^/M (3 '' n,2+ >'4 (2) 

The spin stiffness ps and the spin wave velocity c should 
be viewed as phenomenological parameters to be deter- 
mined either from experiment or from other techniques 
such as spin wave theory or numerical simulations. The 
coupling constant of the model is g — he A/ p s (A is a high 
frequency cutoff). There is a zero temperature quantum 
phase transition at g — g c ~ Att from a phase with long- 
range order (g < <? c , "renormalized classical regime") 
to a phase which exists for g > g c and which is quan- 
tum disordered with only finite spin correlations and no 
static magnetic order. It is now firmly believed that the 



5=1/2 Heisenberg model described by QJ has g < g c . 
Measurements of the correlation length of La 2 Cu04 have 
been fitted extremely well with the QNLerM predictions 
for the renormalized classical regime. 9 

Once holes are added to the Cu0 2 planes, the mag- 
netism becomes rather complex. Fig. ^ summarizes the 
magnetic phase diagram at weak doping concentrations 
of La 2 _ x Sr x Cu0 4 and Y 1 _ x Ca x Ba 2 Cu04ii& Here, we 
concentrate on La 2 _ x Sr x Cu04. For very small Sr con- 
centration, the most dramatic effect is a rapid reduction 
of Tjy with the complete destruction of long-range order 
occurring at a critical doping value of roughly x g ~ 0.02. 
Further, a spin freezing is observed inside the AF phase 
below a temperature Tf which scales linearly with the Sr 
concentration, Tf ~ (815K)x for < x < x g . This spin 
freezing is inferred from a broad distribution of extremely 
slow relaxation times measured with local probes such as 
139 La nuclear quadrupole resonanceii (NQR) and muon 
spin resonance^ (/iSR). Surprisingly, while at higher 
temperatures doping leads to a reduction of the local 
staggered moments, at temperatures lower than about 
30 K the staggered moments recover and at zero tem- 
perature they are practically doping independent and 
approach the value of the undoped compound?ii*i£ see 
the middle panel of Fig. ^ However, the distribution 
of staggered moments is broad at finite doping, with a 
variance which is again simply linear in x, see Fig. ^ 
bottomiifi Both the recovery of the staggered moments 
and the broad distribution of relaxation times are remi- 
niscent of a transverse spin glass state, in which the trans- 
verse spin wave modes of the AF freeze in a static but 
random pattern. These are clear signatures of disorder in 
the weakly doped AF phase. This is further corroborated 
by transport measurements, which show a behavior typi- 
cal for random systems^ At temperatures below ~50 K 
the conductivity roughly follows variable range hopping 
characteristics while at higher temperatures a thermally 
activated conductivity is observed, with activation en- 
ergies of about 19 meVii^ This indicates that the holes 
localize near the randomly distributed Sr donors. 

Both the presence of finite staggered magnetic mo- 
ments and the broad distribution of slow relaxation times 
persist also above x > x g where long-range order is 
destroyed^ Again, there is a recovery of the staggered 
moments at very low temperatures, although the zero 
temperature moment is now slightly smaller than in the 
undoped compound. The x dependence of Tf follows 
now roughly a 1/x scaling. The regime 0.02 < x < 0.05 
is well described as a conventional spin glass (SG) and 
shows characteristic non-ergodic behavior^ The freez- 
ing transition temperature Tf in this regime can thus 
be identified as a SG transition temperature T g . The 
fact that staggered moments persist also above x = 0.02 
is important and excludes the possibility that the tran- 
sition at x = 0.02 is a disordering transition driven by 
quantum fluctuations as described in the QNLerM formu- 
lation above. It is often argued that upon hole doping, 
the reshuffling of the spins by mobile holes leads to en- 
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hanced quantum fluctuations of the spins which would 
eventually drive the spin system past the quantum criti- 
cal point of the QNLer model, driving the AF into a spin 
liquid phase. As the transition at x = 0.02 is not followed 
by a spin liquid phase but rather a SG this scenario does 
not apply for the AF-SG transition. 

Only recently, it was found that the short ranged mag- 
netic order in the SG regime is incommensurate, with 
a maximum of the imaginary part of the susceptibil- 
ity located at the in-plane wave vector (1/2 ± S/y/2, 
1/2 =p S/y/2), in units of 27r/a where a is the Cu lattice 
spacing^£ii£*i£ Here, 5 is the incommensurability which 
roughly follows 5 ~ x. This incommensurability has of- 
ten been interpreted as diagonal stripe formation, even 
though no signatures of a charge modulation were ob- 
served. Rather, all experiments point toward a quenched 
charge distribution and we thus argue that a more likely 
explanation is the formation of short ranged spiral order. 
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FIG. 1: Phase diagram as seen by /i-SR, with data obtained 
from La2- x Sr x Cu04 (open symbols) and Yi_ x Ca x Ba2Cu03 
(closed symbols), p 8 h is the hole concentration, (a) Doping 
dependence of the Neel temperature Tjv, freezing transition 
temperature T/, spin glass transition temperature T g and su- 
perconducting transition temperature T c . (b) Normalized av- 
erage internal field at T=l K. (c) RMS deviation AB at T=l 
K. Fig. from Niedermayer et aim— 

In La2- x Sr x Cu04 static AF moments are strong for 
small x and the holes seem to localize at low temper- 
atures where transport experiments indicate a relatively 
weakly bound hole with a localization length of a few lat- 
tice constants. Thus, one might hope to gain considerable 
insight into these phases by solving the one-hole problem 
first and to proceed from there on. As mentioned in the 
beginning, the understanding of the spin-polaron state 
arising from one hole in an antiferromagnetic background 



is by now quite matureiSii^iiS For the t — J model, the 
bottom of the dressed hole band lies at the zone face cen- 
ters ko = (±7r/2, ±7t/2) and the bandwidth scales with J. 
Because of the presence of two sub-lattices, there exists 
a pseudo-spin degeneracy for each k vector. An impor- 
tant characteristic of the hole wave function is that it 
describes a long-ranged dipolar distortion of the AF or- 
der which arises from a coupling of the spin current car- 
ried by the hole to the magnetization current of the AF 
background^ Relative to the position of the moving hole, 
the Fourier transform of the transverse spin deviations is 
then proportional to (q x +q y )/q 2 ^2, where q = q— (it, it), 
i.e. the staggered magnetization shows a dipolar pattern 
in real space identical to the one produced by an isolated 
ferromagnetic bond, see Sec. Ill Al 

The Sr impurity position, located above the center of 
a Cu plaquette, has a high symmetry and couples to 
both sub-lattices in the same way, so that the pseudo- 
spin degeneracy mentioned above should survive also in 
the bound hole state. The bound hole state is a super- 
position of plane wave states describing the mobile hole. 
For sufficiently weakly bound holes, we expect the main 
weight of the bound hole wave function to remain at wave 
vectors close to ko or equivalent positions, and, depend- 
ing on the relative phases and the weight of these pock- 
ets, dipolar or quadrupolar frustration is associated with 
the localized hole. We note that dipolar frustration was 
also suggested by Aharony et al. for O doped systems, 
caused by a localization of holes in the O site with the 
liberation of one of the spins from the O p 6 state (Si lead- 
ing to an effective ferromagnetic coupling for the two Cu 
spins joint by the O. While the microscopic origin of frus- 
tration in the Aharony model is very different from the 
quantum mechanical one that we assume here, the phe- 
nomenological spin-only model we employ below is not 
sensitive to the microscopic details. In either case, the 
dipole moment of the localized hole state is character- 
ized by two vectors, one in spin- and one in real space. 
The real space vector characterizing the dipole is simply 
the orientation of the ferromagnetic bond in the Aharony 
picture while it is determined by the four coefficients Ck 
and by the equivalent wave vectors of the bound hole 
wave function in the quantum mechanical model. The 
coupling to the spin background is then identical in both 
models. Here we simply assume that the localized hole 
produces dipolar frustration and, rather than relating our 
phenomenological coupling parameters to a microscopic 
model, we derive our parameters from a comparison to 
experiments. As we discuss below, the dipole model can 
quite well explain all the important characteristics of the 
magnetism of the weakly doped AF and SG phase. Let 
us further mention that for Sr doping, it was proposed 
that a chiral spin current is induced on the four Cu sites 
closest to the Sr impurity which leads to a Skyrmion-likc 
distortion of the AF, where the mechanism of frustra- 
tion is again the coupling between spin and background 
magnetization currents^ 

In section [H] we introduce the dipolar frustration 
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model, summarize the main results of previous studies 
on this model, and discuss how they compare with ex- 
periments. In section Hill we first derive an extension of 
the model to allow for non-collinear correlations which 
arise from dipole ordering. We perform a RG calcula- 
tion to understand the influence of disorder and discuss 
the importance of topological defects of the spin texture. 
Finally, in section Hvl our results are compared with neu- 
tron scattering data on the SG phase of La 2 - x Sr x Cu04. 
We find that the SG phase can be described as a strongly 
disordered spiral phase in which topological defects pro- 
liferate. 



II. THE AF PHASE AND DIPOLAR 
FRUSTRATION MODELS 

We briefly sketch here the basis of the dipolar frus- 
tration model and the results of previous studies of this 
model in the collinear limit. The model as presented in 
this section is applicable only for the antifcrromagnetic 
phase in which the dipoles do not have a preferred di- 
rection. At high temperatures, the collinear theory can 
then be used. We will show in the next section however, 
that the collinear model is not able to describe the low 
temperature and/or strong disorder regime, where non- 
collinear behavior emerges. 

In the dipole model, it is assumed that each localized 
hole produces dipolar frustration. It is then possible to 
study the magnetism of the hole doped materials com- 
pletely ignoring the charge degrees of freedom and to 
work with the spin sector only. Further, as there are clear 
indications of static AF correlations for x < 0.05, the an- 
tiferromagnet should be well described within the renor- 
malized classical regime of the QNLctM. In this regime, 
quantum fluctuations simply lead to a renormalization of 
the coupling constant of the classical model. A classical 
model should thus suffice to describe the relevant physics 
in the AF and SG regime. 



Ferromagnetic bonds as an example of dipolar 
frustration 



Dipolar frustration was first discussed in the general 
context of insulating spin glasses by Villain^ The sim- 
plest way of producing dipolar spin textures is by placing 
a ferromagnetic bond in an otherwise AF magnet, whose 
order parameter we denote by n. At a distance x away 
from the ferromagnetic bond, this leads to a deviation of 
the Neel order Sn ~ f^x^/x 2 . Here, f M is a vector both 
in spin and lattice space, where fi — 1,2 are the indices 
of the 2D lattice vector. The spin part corresponds to 
the local ferromagnetic moment (with _L n) produced 
by the bond and the lattice part corresponds to the ori- 
entation of the bond on the lattice (see Fig- EJ - This can 
be easily derived in a harmonic continuum approxima- 
tion, where the energy density of the magnet away from 



the impurity is proportional to [<9 p ((5n)] 2 and the classi- 
cal equation of motion is V 2 (<5n) = 0. For any impurity 
distribution, the solution for i5n can thus be written in 
a multipolc expansion. As the monopole moment is en- 
ergetically too costly 2 ^ the lowest order contributions, 
consistent with the symmetry of the one-bond problem, 
are dipolar. 
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Dipolar distortions produced by a ferromagnetic 



B. Collinear Model 

Because of the long-range nature of dipolar frustration, 
a continuum field theory, such as a (classical) non-linear 
a- model (NLcrM), should be well suited for a treatment of 
this problem. While the dipole spin structure discussed 
above is a solution of the harmonic theory, it is not a 
solution of the 2D NLcrM. Nonetheless one can study 
the dipole model within the NLerM, if one introduces the 
dipolar frustration through a minimal coupling scheme. 
As mentioned in Ref. [24], the dipolar frustration can 
be reproduced (on the harmonic level) via a coupling of 
the dipoles to the gradient of the order parameter n of 
the NLcrM. Thus, within a NLcrM approach, the reduced 
Hamiltonian of the model can be written as^&Si (the fac- 
tor f3 = T _1 is included in the Hamiltonian and we set 
k B =l) 



H co l — 



Ps_ 

2T 



J d 2 x(d M nf 



T 



d 2 x f„ • <9„,n 



(3) 



where n 2 = 1, p s is the spin stiffness (renormalized by 
quantum fluctuations), T the temperature, n a three 
component unit vector representing the local staggered 
moment and f M is a field representing the dipoles. We did 
not include here small corrections which lower the spin 
symmetry from Heisenberg to XY or Ising. While these 
are known to be present both in the undoped and weakly 
doped compounds^, they have a very small character- 
istic energy scale and, as a first approximation, we set 
them to zero. Note however that these terms dominate 
the static magnetic susceptibility near the Neel transi- 
tion. For a random distribution of localized dipoles we 
write 



x, a, 



(xi)Mi 



(4) 
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where the sum is over the impurity sites, a, are lat- 
tice unit vectors, Mj unit vectors in spin space, and M. 
measures the strength of the dipoles. While there is no 
dipole-dipole interaction term in Eq. ©, fluctuations of 
the n field generate a spin wave mediated interaction. 
This can be seen once short scale fluctuations are inte- 
grated out under a renormalization procedure^ An in- 
tegration over the short scale fluctuations up to a scale 
L> 1/ yfx (but L <§; £ where £ is the 2D spin correlation 
length) leads to an effective interaction term of the form 



mm}} 



PsM 2 
2T 



J ^ Mi 



Mi 



with 



Jin _ I Z n cxi • ci 



x 2 



> 



(5) 



(6) 



and Xij — Xj — Xj. Thus, for an average separation of 
dipoles ~ there is a random interaction among 

dipoles with a characteristic energy U ~ p s 7M 2 x/47r. ft 
was further shownS^ that at high temperatures, U <T, 
the presence of dipoles lead to a renormalized effective 
stiffness p e s = p s (l — U/T). Thus, the correlation length 
at high temperatures (and small x) has the form 



exp 



2irp s 2wp s U 



T 



rp 2 



(7) 



This result agrees to lowest order in x with that ob- 
tained by Cherepanov et alw^ in a related renormaliza- 
tion group (RG) calculation where they calculated p c s 
up to order x 2 . ^From a comparison with correlation 
lengths obtained from neutron scattering data at high 
temperatures, they estimated U ~ 20,052;. The doping 
dependence of Tn was also found to be consistent with 
the dipole modeli^i 
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FIG. 3: Ri = (Ti) -1 data from 139 La NQR relaxation mea- 
surements for La2- x Sr x Cu04 and various x < 0.02, from Ref. 

A second independent test of the value of U is to 
consider the spin relaxation times inside the AF phase. 



This can be understood already within the theoretical 
framework just presented, using arguments similar to 
those from Ref3i where spin relaxation has been dis- 
cussed within a slightly different frustration model. The 
relaxation rates inside the AF phase can be explained 
within the dipole theory if one assumes that the relax- 
ation is driven by the interaction among dipoles and 
hence controlled by the parameter U. At temperatures 
well above the actual freezing temperature, an Arrhenius 
law is observed, with a characteristic energy E = 8.91/ ~ 
7250Kir, see Fig.|3 The above estimate of U correctly re- 
produces the linear scaling of the relaxation energy with 
x and also gives a good estimate for the slope. With 
U = 20p s x, p s - 24 meV 24 one obtains U - 5500Kx. 
Considering that this is a very rough approximation, the 
value is not too far off from the experimental one. We 
mention further that the linear scaling of the width of the 
distribution of local staggered moments is also consistent 
with a dipole modeliSi 



III. NON-COLLINEAR CORRELATIONS AND 
DIPOLE ORDERING 

While the dipole model presented above can well ex- 
plain the temperature and doping dependence of the cor- 
relation length not just in the AF but also, to some ex- 
tent, in the SG regime^i theoretical investigations of the 
model have always predicted (or rather assumed) short 
ranged commensurate antiferromagnetism. The recent 
observation of incommensurate (1C) correlations for the 
regime 0.02 < x < 0.05 requires therefore a new approach 
to the SG phased 

As a possible explanation for the presence of IC cor- 
relations, a disordered striped phase has been proposed, 
similar to the ordered striped phase found near x ~ 1/8. 
While there is indeed an instability in the striped phase 
toward a disordered phase at low x^ it is unlikely that 
the stripes will survive in presence of strong disorder. In 
fact, recent numerical simulations of Shraimann-Siggia 
dipoles with disorder have shown that the latter leads to 
a destruction of the stripe phased 

In the spin glass regime, there are two competing 
length scales. The first is related to the average sepa- 
ration between disorder centers (Sr ions) £d which scales 
as id ~ 1/ y/x. The other is the scaling of the periodicity 
i s associated with the incommensurability, which scales 
as £ s ~ I jx. For small x, id <§; i s - In a stripe scenario the 
charge distribution would also have a periodicity which 
scales with £ s . Thus, in a striped phase the charge can 
not take full advantage of the disorder. The stripes must 
either break up into short segments or reduce their on- 
stripe charge density considerably to take advantage of 
the disorder potential. Instead we propose here a theory 
in which the charges are completely disordered and the 
incommensurability exists only in the spin sector. Then, 
there is no conflict between the two scales i s and id as 
i s relates only to the spins whereas id characterizes the 
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charge distribution. 

Notice that even in the case that short segments of 
stripes should be present, these stripes would lose their 
anti-phase domain wall character and instead act like a 
row of ferromagnetic bonds, again causing dipolar frus- 
tration. Thus, the theory we present here applies both 
to the case of localized hole states which produce dipo- 
lar frustration as it does to a system of randomly placed 
stripe segments. We view the scenario of localized holes 
however as the more plausible one. 




FIG. 4: The order parameter of collinear magnets, which are 
invariant under rotations around the collinear axis, can be 
represented by a unit vector (left), whereas non-collinear or- 
der parameters require three orthonormal vectors (right). 



A. Dipole ordering 

It is easy to sec how the dipole model can lead to IC 
correlations^ The Hamiltonian Eq. favors the for- 
mation of a spiral phase, with a non-zero average twist 
d^n of the AF order and a simultaneous alignment of the 
dipoles, (f M ) 7^ 0, as long as the lattice and spin degrees of 
freedom of dipoles are annealed and free to orient them- 
selves. The lattice position of the Sr dopants (located 
above the center of a Cu plaquette), which pin the holes, 
suggests that this freedom indeed exists. We emphasize 
that a spatially homogeneous distribution of dipoles is 
not required for the formation of spiral correlations. 

The preferred orientation of the lattice part of the f M 
vector is determined by the nature of the localized hole 
state and therefore should reflect the symmetries of the 
underlying lattice. Thus a discrete set of favored lattice 
vectors for the formation of the spiral exists. The a-b (or 
square lattice) symmetry breaking associated with the 
formation of spiral correlations can therefore have truly 
long range order. The continuous symmetry of spin space 
on the other hand inhibits long-range magnetic order in 
the 2D system for either finite temperatures or disor- 
der. The experimental observation of a macroscopic a-b 
asymmetry^ but very short spin correlation lengths thus 
clearly motivate the study of the dipole model. 



The spin and lattice symmetries are decoupled and inde- 
pendent for the collinear AF. A spiral ground state on the 
other hand breaks the 0(3) spin symmetry completely. 
Moreover, in a spiral state the lattice symmetries and the 
spin symmetries arc no longer decoupled and the order 
parameter space of such a state becomes more involved. 

For spirals, the combined symmetry of lattice and the 
spin space is 0(3)xO(2). As discussed in detail by Azaria 
et al.jli> the coupling of the spin and lattice degrees of 
freedom in frustrated spin systems leads to an order pa- 
rameter which results from a symmetry breaking of the 
combined lattice and spin degrees of freedom and is in 
general of the form 0(3) x 0(q)/0(q) where q depends 
on the symmetries of the lattice. For a spiral phase, one 
finds2£ q = 2. 

A convenient representation of the order parameter is 
in terms of orthonormal n^, k — 1...3, with n^,n a q = 
Sk q - Klee and Muramatsu^ have derived a continuum 
field theory for the rife order parameters from the lattice 
Heisenberg model Eq. QJ, assuming an IC spiral state 
with an ordering wave vector kg = (-,—)+ qg. Here, 
qs measures the deviation from the commensurate AF 
wave vector, see Fig. |SJ The lattice spins Si at sites 
can be parametrized in a spiral configuration with the 
use of the n£ as (with 113 = ni x 112) 

Si/S = ri! cos(k s • r 2 ) - n 2 sin(k s • r t ). (8) 



B. Continuum description of spiral phases 

We here investigate the dipole model allowing for the 
presence of non-zero ordered moments but assume a ran- 
dom spatial distribution of the dipoles. First, however, 
we need a proper theoretical description of the homoge- 
neous spiral phase. 

In collinear magnets, the rotational 0(3) symmetry 
of the system is broken down to a ground state with 
0(2) symmetry, as rotations around the magnetization 
axis leave the ground state invariant (this is schemati- 
cally shown on the left hand side of Fig. |3J. The or- 
der parameter of collinear magnets is then an element 
of 0(3)/0(2). This group is isomorphic to the group 
of three dimensional unit vectors n, which is the repre- 
sentation used in the Hamiltonian Eq. J3J. Further, in 
absence of dipoles, the Hamiltonian Eq. ((3J) is invariant 
with respect to 0(2) rotations of the lattice variables. 



^\ I x\ f v\ I , 

x\ f v\ \ f , ^\ t A 

A f — v\ I v\ f 
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FIG. 5: Spin texture of an AF spiral. 

A perfectly ordered spiral is described by Eq. © with 
constant, i.e. space independent, n^. To allow for spatial 
fluctuations of the spins around the spiral order, Klee and 
Muramatsu introduced a slowly varying field L via£2iiS 



S(r t 
S 



N + aL 



VI + 2aN • L + a 2 L 2 
= N + a [L - (N • L)] - ■ 



(N-L)L + -L 2 N 
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^(N-L) 2 N 



+ 0(a 3 ), 



(9) 



where N = ni cos(ks -r,-) — 112 sin(ks -r^) with now slowly 
fluctuating fields n^. The continuum theory can then be 
found upon expressing in the lattice Heisenberg model 
the spin operators in terms of the and L fields, ex- 
panding the terms up to order a 2 and taking the limit 
a — > in the end. After integrating out the L fields, one 
finds an effective Hamiltonian which can be written in the 
classical limit in the general forrn^ 2 . (again we include the 
factor {3 — T^ 1 into H) 

H=\jd*x Pkiiid^k) 2 + s» J d 2 x ri! • ^n 2 . (10) 

This description is valid for length scales larger than 
Iqsl^ 1 . The stiffnesses of the order parameter are 
given initially by p lfl = p 2f _ l = JS 2 cos(q Sll a)/(2T) and 
P3 M ~ 0, but will change under a renormalization of 
the model. We will ignore for the most part the small 
anisotropy (of order q^a 2 ) in the stiffnesses php, and 
just write pk- The vector s is to lowest order given by 
s = Jqg/T. The term with the s M pre-factor makes 
this Hamiltonian unstable, which simply expresses the 
fact that the pure Heisenberg model does not support 
a spiral phase ground state. The term will however 
be canceled by a similar term originating from the cou- 
pling of the spins to the ordered fraction of the dipoles, 
relating the incommensurability self-consistently to the 
ordered moment of the dipoles. In other words, the or- 
dered dipoles stabilize the spiral phase, as expected. 

It must be stressed that because the continuum model 
is only valid at length scales larger than the period of the 
IC structure, there is a relatively large uncertainty in the 
estimates of the Php,- There is always a fundamental prob- 
lem in relating the continuum model parameters to those 
of the original microscopic lattice model, but in this case 
this problem is especially severe. The continuum model 
parameters must be obtained from an average over one 
period of the spiral which, for small incommensurabili- 
ties, can be rather large. Thus, the above estimates for 
the pfc^'s should be taken with care. 



C. Disorder coupling: a gauge glass model 

We now must include the coupling of the dipolar frus- 
tration centers to the spiral order parameter. While there 
is no microscopic derivation of this coupling at hand, the 
fact that the coupling in the collinear model can be ex- 
pressed within a minimum coupling scheme allows for a 
simple generalization of the model to non-collinear spin 
states. We first observe that the ordering wave vector 
of the spiral, qs, is entirely determined by the aver- 
age orientation of the dipoles. Similarly, local variations 
of the density or orientation of the dipoles should also 
modify the local ordering wave vector. Further, to re- 
produce the strong canting produced by the dipoles, the 



coupling should be of first order in the spatial derivative 
of the spiral order parameter. To generate the frustra- 
tion produced by the dipoles we thus introduce a minimal 
coupling 34 in the first term of Eq. (|10fl . i.e. we replace 
(c^rifc) 2 with [(<9 M — zB M • L)rife] 2 where B M is a random 
gauge field, representing the dipoles. The components of 
L are 3x3 matrix representations of angular momenta 
which generate rotations about the three spin axes, with 



iBa • L rife = B M x n fc . 



(11) 



This coupling has the advantage of relative simplicity 
combined with a clear physical interpretation: the dipo- 
lar fields define the locally preferred wave vector of the 
spiral, and fluctuations of the dipole fields lead to fluc- 
tuations of the wave vector. Further, it reproduces the 
correct form of the dipole coupling in the collinear limit, 



as shown below. Let us write B M = [B^J^, + Q u so that 



[Q/Jd = u > wnere [■ ■ -Id i s t ne disorder average. We 
then obtain the following Hamiltonian for the spiral in 
presence of disorder, 



i J d 2 xp ktl (d ll n k ) 2 + J d 2 xp fe a M n fc -Q M xn fe , (12) 



H = - 
2 



where the ordered part of the dipole field cancels the 
second term in Eq. (|10|) . Thus, 



Pk^dftiik ■ [B M ] D x rife + s M ni ■ 9 M n 2 = 0. 



(13) 



As qs oc s, this equation relates the incommensurability 
linearly to the density of ordered dipoles. The remain- 
ing part of the dipole field, Q M , is a quenched variable 
with zero mean and we assume Gaussian short ranged 
statistics, 



[Q£(x)Q»(y)] =A5(x-y) S ab 6^. 



(14) 



In absence of disorder, the Hamiltonian defined by 
Eq. (H2) has the desired (9(3) x (3(2)/0(2) symmetry. 
The 0(3) symmetry is associated with the spin indices a 
of the n k , while the 0(2) symmetry is associated with the 
lattice indices k and arises because p\^ = p2^- Hence, the 
equality pi^ = p 2 ^ is directly related and enforced by the 
symmetries of the spiral. Note that if all p klx are identi- 
cal, the lattice symmetry is enhanced to 0(3). We further 
see now, that the model reduces to the collinear model 
Eq. m the case p\^, = with p$ = p s /T, 113 = n and 
f M = Q M x n. Unfortunately it is not possible to reach 
the collinear limit by sending qs — > 0. The reason is that 
the parameters p kjl are, within the approximation used 
in their derivation, independent of the size of the unit cell 
of the spiral, i.e. in the limit qs —* 0, the unit cell size 
diverges while the parameters pfe M remain unaffected. 

The model defined by Eq. I|12l) is in fact far more gen- 
eral than its derivation might suggest. In absence of dis- 
order it is applicable to other types of frustrated spin 
systems with a non-collinear ground-state, such as e. g. 
the Heisenberg model on a triangular latticei^iiS 3 - 3 ^ It 
is conceivable, that certain types of randomness in such 



8 



lattices may be well described by the disorder coupling 
employed here. More importantly, the model Eq. i|12|) 
can be viewed as a general model to investigate diluted 
spin glasses, in which a spin system is frustrated by a 
small number of impurities. There have been investiga- 
tions of similar models of spin glasses in the past, most 
notably by Hertz^4 which however did not account for 
non-collinear correlations which are known to be essen- 
tial in spin glasses.— Our approach has the appeal that it 
can interpolate between collinear and non-collinear states 
and thus offers the possibility to study the transition from 
an ordered collinear magnet to a disordered non-collinear 
one continuously. 

D. Renormalization 

We now investigate the renormalization of the model 
under a change of scale, with the objective to understand 
the influence of the dipoles on the correlation length of 
the model. For carrying out the RG calculation, it is of 
advantage to use a SU{2) representation of the model2& 
(see also App. lA"ll. We therefore write 

1 

-1 

2 



= . tr KV fe 9 1 



(15) 



where a a are Pauli matrices and g 6 SU(2). We further 
introduce the fields^ 



4 = ^-tr [a a g-%g] , 



(16) 



which are related to the first spatial derivatives of n/. 
through dfjjil = 2t i ^ k A l il n a y Eq. (|12fl then acquires the 
form, 



H 



bA 7 : 



4n>£<#. (17) 



where t^ 1 = 2(p lf _ l +p 3fi ) and b = (p lp -p 3 ^)/(pip +P3/0- 
At the point b = the symmetry is enhanced to 0(3) x 
0(3)/0(3) ~ 0(4)/0(3) while at b = -1 the model is 
collinear. For spirals, we have initially 6 = 1. 

We first discuss the dimensional scaling behavior of 
the model ((121 I17|l . We assign the dimension —1 to each 
spatial dimension so <9 M has dimension 1. It follows that 
the A M fields have a scaling dimension of 1. The scaling 
dimension of the first term in Eqs. I]12l 117(1 is then 2 — d 
where here d = 2. Thus, this term is marginal and an 
RG analysis is required to study the scaling of the t tll b 
parameters. Local terms containing more than two A M 
terms have positive dimensions and are irrelevant. Hence, 
such terms, while they are generated in the perturbative 
expansion we discuss below, need not be considered. 

As was pointed out in Ref . [24] , for the disorder choice 
(|14fl the model defined by Eq. (J3J) has a lower critical di- 
mension of two and is thus renormalizable in two dimen- 
sions, as can be shown with a general Imry-Ma type ar- 
gument. The same argument can be used for the present 



model. The disorder coupling in Eq. 1(12(1 can be rewrit- 
ten in momentum space as a random field coupling of the 
form 



/ 



d 2 q 



njfe(-q) • hfe(q); 



(2tt)2 

hfc(q) = ipkv% J d2yi (Qm x n fc)e 4q ' x 

where the random fields (q) have disorder correlations 

^(q)^;(qOl„cx%- 



D 



with a momentum dependence 

q')|q| e with 8 = 2. According to general arguments by 
Imry and MajS in models with continuous symmetries 
random fields will destroy long-range order as long as 
d < 4—0. This implies that in our case d = 2 is the 
lower critical dimension^ and a renormalization group 
analysis of both the stiffness and the disorder coupling is 
required. 

We now derive the one-loop RG equations. For this, we 
split the original SU (2) field g into slow and fast modes, 
g = gexp(i ip a a a ) and trace out the fast modes ip a which 
have fluctuations in the range [A -1 , 1], where we set the 
original UV cutoff equal to 1. For the one-loop calcu- 
lation, we need an expansion of Eq. (|17l) up to second 
order in ip a (higher order terms will only contribute at 
higher loop order of the RG). For the fields and A M 
the expansion reads (see App. for more details) 

4=4 + + e ijkV Pd^ k + 2e^4 - 24 ^ 
+ 2A M -<? <p l +0{</), 

where 



K jk 2 1 V _ 2 2 ' 



The expansion of the energy functional 1(17(1 reads 



H = 
with 

H 



d 2 x 



6 lA! 



H, 



rO 



H^ + Hp (18) 



d 2 xpk^ijk e a bc A^hlQ^, 
Hi + H c i + H C 2 



p — H i + Hi + H3 



H, 



c3 



H c4 . 



The first two terms in the expansion of H have exactly 
the same form as the original functional l|17l) . but are 
now functionals of the slow fields. H v is quadratic in ip 
and has the form 



H uy = — I d 2 x 



H\ . . .H\ are generated by the first term in Eq. (fTTjl and 

are given by 

H x = 2t- x [ d 2 x A i li d^i P k e ijk (1 - b5 iz + 2b5 jz ) , 
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H 2 = 21- 1 J d 2 x d^Ai (1 + b5 lz ) , 
H 3 = ibt- 1 J cPxe^AtfA*, 
Hi = ibt^ 1 J d 2 x 

The coupling term in Eq. Ijl7|l produces the H c i . 
terms, 

H cl = Ajd^ Pk ^ k e abc 



■ H, 



cA 



^kimnjh^A 1 ^ 



e jimn%nf n A l fi + tii m n a jh c k A™ 



H c2 = 2 J d 2 x p k ^e ijk € abc d^ip' l h^h L k Q^ 
H c3 = 2jd^ Pk ,e vk e abc 



4 (^m + n c k Rt) 



+ 2h«h% (A^ 

+ a I \i ~ a ~ c I a m ~ c ~ a 



A^e Um e jm h a q n'Cj <^V 



H c i — 2 J d 2 x p k ^ij k e abc [2 (e k i m n a J n c m 

The integration over the fast tp fields is performed with 
J exp(-F v ) exp(-F p ) = ^ J exp(-H v ) 

where T is obtained from a cumulant expansion 

_ T = ln /Pb i ]exp(-g y )exp(-g p ) 
/2%*]exp(-#„) 

71=1 

and (. . .) indicates that only connected diagra 
to be considered. 



rams are 



With ^ = lnA and 



C x (0) 



Air 



■In A 



one then finds the RG equations 



d_l_ 

d£ t s 
d_b_ 

dIT s 



l-b _ {2-b + b 2 )X 

2tt ~ Ant s 
(3 + 6)6 (5 + 6)6A 



2tt 



Ant s 



This yields 
dl ts 



l-b 



2tt 



-ti + 



2-b + b 2 



-in 



At., 



6(1 + 6) 6(l + 6)(3-6) 



47T 



A. 



(20) 
(21) 



For A = 0, these equations describe the RG of a clean 
spiral^ while for the collinear point b = — 1, the equa- 
tions reproduce the RG of the stiffness for disordered 
collinear models^ From Eq. (|21() it is seen that there 
are two fixed points for b (the asymptotic freedom of the 
model prevents a true fixed point in 2D as t s always di- 
verges). The collinear point b = — 1 is unstable whereas 
6 = is stable, irrespective of the disorder. The RG flow 
of t s and b is shown in Fig. [fj] for A = 0. The flow does 
not change qualitatively for finite A as long as A -C t s . 
Hence, the coupling to weak disorder does not lead to 
any new fixed points, although the disorder renormalizcs 
the stiffness. 
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1. Renormalization of the spin stiffness 

We ignore the (small) anisotropy of the parameter 
and simply use the isotropic mean t s = \Jt\ti in the RG 
analysis below. We collect all terms in the perturbative 
expansion which are bilinear in A 1 ^. After performing the 
disorder average of T , the renormalized stiffnesses of the 
A 1 ^ fields is found to be (see Add. iGland ID 1|) 



b=-1 



2(1-6) (2-6 + 6 2 )A 



+ 



26(3 + 6) 6(5 + 6)A 



C x (0), 



t 2 



C x {0). 



b=0 
b 



b=l 



FIG. 6: RG flow of t s and b for A = 0. For any b > -1, the 
flow is toward 6 = 0. 



2. Renormalization of disorder coupling 

As we discuss below, the renormalization of A is given 
by terms proportional to Xt s and A 2 . As the disorder 
enters the renormalization of t s only in the combination 
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Xt s (see Eq. JUJl), we can neglect the renormalization 
of A altogether for t s ^> A, i. e. at high temperatures 
(we have t s oc T/J). However, for low temperatures the 
renormalization of A must be taken into account. To cal- 
culate the renormalization of the disorder we follow the 
approach used in Ref. [24]. In this approach, the renor- 
malized disorder variance is defined by the variance of 
all terms in the perturbative expansion which couple to 
the quenched disorder fields and are linear in A^. Note 
however that there exists no symmetry argument which 
guarantees that the functional form of the disorder cou- 
pling remains unchanged under the RG. It is thus pos- 
sible that new disorder terms are generated so that a 
simple renormalization of A is not sufficient. This is in- 
deed the situation we encounter for general 6 7^ and 
discuss in more detail below, where we find the gener- 
ation of new coupling terms at order A 2 . To find the 
complete renormalization of the model one would have to 
include all generated new terms into the original model, 
which is a rather laborious process which we did not pur- 
sue. Nonetheless, as we have just shown above, there are 
only two possible fixed points even in absence of disorder, 
6 = and 6 = — 1. Rather than trying to categorize all 
possible disorder couplings, we therefore focus on a dis- 
cussion of the RG of the disorder near these two possible 
fixed points and discuss their stability under the flow. 

We begin with the collinear case, 6 = — 1. In this case, 
the renormalized variance of the terms linear in A 1 ^ is 
given by (see App.EHand App.El Eq. (JDT3)) 



¥ 1 d2x 
1 



1 t s In A 

7T 



- — A In A 

2tt 



M ) 2 J-AlnA 

2tt 



(22) 



What is evident from this result is that the renormal- 
ized disorder coupling is no longer of the original form 
Pkd^rik ■ Q M x rife. Such a coupling has a variance which 



includes a prefactor of (1 + 6) of [A^J ■ According to 

Eq. I|21l) . 6 = — 1 is not changed under the influence of 
the original disorder coupling. A renormalization which 
retains the form of the original coupling can then not lead 
to a renormalized disorder variance with a finite prefac- 
tor of a t b = — 1. Such a term is however present 

in Eq. I|22ll we conclude that a new type of disorder cou- 
pling is generated at 6 = —1. This is perhaps easier 
to see in Fourier space, where the original disorder cou- 
pling can be written as a correlated random field coupling 
rifc(— q) ■ hfc(q), see Eq. (|18|) . For the original minimal 
coupling one has h/-(q) cx p^ and thus, in the collinear 
limit 6 = — 1 (or pi = 0) , only 113 is affected by this cou- 
pling. We can then interpret the finite prefactor of the 
(j4*J term in the disorder variance as the generation of 
correlated fields which couple also to ni j2 even at 6 = —1. 
It is evident that such a coupling will drive the system 
away from 6 = — 1 and thus destroy the collinear fixed 
point. Thus, even if the original AF order is collinear 



(i.e. in absence of dipole ordering), the disorder drives 
the system to a non-collinear state. An analysis which 
pre-supposes collinear order is thus not valid in the pres- 
ence of dipoles and cannot describe the low temperature 
regime correctly. Physically, one would also expect the 
appearance of non-collinearity. The random canting of 
spins leads to a random local deviation of the spins from 
the ordering axis and thus destroys the remaining 0(2) 
spin symmetry of the collinear model. 

To make contact with the RG result obtained from the 
collinear model in Ref. [24], we note that we can repro- 
duce the result Cherepanov et al. obtained for the dis- 
order renormalization if we ignore non-collinear modes. 
We can then define the renormalization of A just by the 
terms which are present in a purely collinear theory, i.e. 
2 / _ \ 21 



by the 



term in Eq. Then 



a a 



2A A 2 
TrtT ~ 2-KtV 



which, using Eq. gDJ leads to 

9 A - 3 A 2 
d£ X ~ 2^ A • 



(23) 



(24) 



This, together with Eq. (|20f) are the RG equations found 
in Ref. [24] (note that our stiffness t s differs from the 
stiffness t used in Ref. [24] by a factor two) . We emphasize 
that this result ignores the role of non-collinearity in the 
problem. 

We now turn to the point 6 = 0, the only remaining 
possible fixed point of the model. At this highest symme- 
try point we find that no new coupling terms are gener- 
ated. The variance of the renormalized disorder coupling 
takes the form 



4t s + 3A 

47T 



In A 



Thus, 



a a 



1 A 3 A 2 

IT t* 47T tl 



(25) 



(26) 



which yields the RG equation, valid for 6 = but any 
initial ratio of X/t s , 



— A - — 

dl 4tt' 



(27) 



Using Eq. I|20|) . we can simplify this through z = f s + A/2 
to get 



d_ 



1 



(28) 



So for 6 = the presence of disorder leads to an ad- 
ditive renormalization of the stiffness, t s — > t s + A/2. 
In presence of any amount of disorder, the IC correla- 
tion length £ at T — is finite, as can be inferred from 
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an integration of the RG equation with 6 = 0, yielding 
£ oc exp(C(i s0 + A /2)~ L ) with some cutoff dependent 
constant C. Thus, even at T = 0, £ oc exp(2CAg 1 ) 
is finite. While the disorder scales to strong coupling, 
the relative disorder strength with respect to the stiff- 
ness, X/t s , always scales to zero so that at long wave- 
lengths the disorder becomes less relevant. This is sur- 
prisingly different to the situation with b = — 1 fixed)24 
where the ratio X/t s was found to diverge below a cer- 
tain initial value of Xo/t s0 which was interpreted as the 
scaling toward a new disorder dominated regime. Thus, 
if one correctly takes into account the non-collinearity, 
this disorder dominated phase disappears. The absence 
of a sharp cross over from a weak disorder to a strong 
disorder regime is certainly surprising, especially as the 
experiments clearly observe a transition into a spin glass 
phase at a finite temperature^ The finite temperature 
transition may be related to the presence of inter-layer 
coupling. We argue below, however, that topological de- 
fects can alter the RG behavior considerably and may 
be a more natural explanation for the appearance of a 
strong disorder regime. 



E. Topological defects: saddle point treatment 

The RG results presented above do not take into ac- 
count topological defects^ of the spiral as only spin 
waves excitations enter the calculation. As is well known 
from XY spin models, topological defects can play an im- 
portant role and drive finite temperature transitions. 40 
The neglect of topological defects has been a source of 
criticism toward the NLctM approach to frustrated mag- 
nets, which gives controversial results for e = 1, 2 in an e 
expansion around D = 2+e dimensions 4i For two dimen- 
sional systems, the NLctM results were however found to 
be in very good agreement with numerical simulations as 
long as the temperatures were sufficiently low 4^ Only at 
higher temperatures, a deviation from the NLctM predic- 
tions for the temperature dependence of the correlation 
length was observed which was attributed to the appear- 
ance of isolated topological defects. In the numerical sim- 
ulations the high temperature region showed some resem- 
blance to the high temperature region of XY-models 4 ^ 
which indicates that this region is characterized by free 
defects. However, at present a good understanding of the 
influence of such defects in non-collinear systems is still 
lacking 

The topological defects of spirals have their origin in 
the chiral degeneracy of the spiral, i.e. the spiral can 
turn clock- or anti-clock wise^i At a topological defect, 
the spiral changes its chirality. As the chirality takes only 
two possible values, the defects are Zi defects. 

It is then straightforward to find topological defect so- 
lutions of the saddle point equations of a clean spiral^ 
The saddle point equations can be obtained from the per- 
turbative expansion of the energy density, Eq. I|18ll9fl . 
One finds that extremal solutions must satisfy for each 



j = x,y,z the equations 

(l + b5 jz )d fl Ai = 2be zjk AlA* 



(29) 



where j is not summed over. For b > — 1 one finds solu- 
tions of the formSi 



(30) 



g s (x) = cxp ( ~mV°*(x) 



where m is a space independent unit vector and 4"(x) 
a scalar function. With this Ansatz, one has A*(x.) = 



^to'c^^x) and thus, upon insertion into Eq. I|29|) . one 
finds for m and W the equations (j is again not summed 
over) 

(1 + b6 jz ) m^tf(x) = be zjk m z m k (d M *(x)) 2 . (31) 

The weight of the configuration described by g s is given 
by (we set t^=t s ) 



H[g s 



1 



b(A z 



l + b(m z y / cfx (d^y. (32) 



We see that for b < 0, the energy is minimized for 
(to 2 ) 2 = 1 whereas for b > the vector m is preferably 
orientated within the x-y plane with m z = 0. For both 
cases, Eqs. Ij31|l reduce to the two dimensional Laplace 
equation V 2x I'(x) = 0. This equation allows for topolog- 
ical defect solutions with ^(x,y) = arctan(j//x). In the 
top of Figs. □ and the spin distribution around isolated 
defects is shown for both b < and b > 0. Using Eq. I|32|) 
one finds that the energy of a topological defect solution 
^(x,y) diverges logarithmically with the linear system 
size i?, 



1 + (m z fb , 
BE = — — '— Trlni?. 

2t, 



(33) 



Because of this logarithmic divergence of the energy, 
isolated defects are not present in absence of disorder 
and at sufficiently low temperatures. It can also be 
shown,— that a bound state of defect pairs, described by 



9 = g s i9s2 with g sli2 = exp 



;mi,2 • <r arctan 



/ y-yi,: 

\ x — xi 7 



has a finite energy if mi + m2 = 0. Therefore, while iso- 
lated defects may be absent, defect pairs will be present 
at any finite temperature. Figs. [7| and [H] (bottom) show 
such a pair of topological defects for b < and b > 0, 
respectively. 

This situation is reminiscent of the one encountered in 
the XY model where at low temperatures also only de- 
fect pairs are present. The pairs unbind at the critical 
Kosterlitz-Thouless temperature. An unbinding of de- 
fects at a critical temperature or critical disorder strength 
is also expected in the present model. The topological 
defects of the spiral differ however in important aspects 
from those of the XY model. Spiral defects have a Z2 
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FIG. 7: Single topological defect (top) and topological defect 
pair (bottom) of a spiral with b < (small scale AF fluctua- 
tions are not shown). 



FIG. 8: Single topological defect (top) and topological defect 
pair (bottom) of a spiral with b > (small scale AF fluctua- 
tions are not shown). 



charge while XY defects have Z charges. More impor- 
tantly, as the present model possesses asymptotic free- 
dom, it has a finite correlation length £ at any finite 
temperature even in absence of free defects. This im- 
plies that the logarithmic divergence in Eq. (|33fl appears 
only up to a scale R < £. It is therefore not clear how 
a defect-unbinding would affect the system. A transition 
from a phase with algebraically decaying spin correlations 
to a phase which shows an exponential decay, as occurs 
in XY models, is clearly ruled out. While in XY models 
topological defects can be relatively easily incorporated 
into the analysis because they can be decoupled from the 
spin waves, this is not the case for frustrated Heisenberg 
models. If fluctuations around the saddle point solution 
are taken into account, the defects of spirals couple to the 
spin waves already at second order in an expansion in the 
fluctuations^. These difficulties have to date prevented 
a good understanding of defect unbinding in frustrated 
systems. 

A comparison to XY models is nonetheless quite illu- 



minating. The kind of disorder coupling we have used for 
the spiral phase is closely related in spirit to XY mod- 
els with randomly fluctuating phases, where the disorder 
is also introduced in the form of a fluctuating gauge. 44 
If one ignores vortices, the influence of the disorder was 
shown to amount to a simple renormalization of the spin 
stiffness, at all orders in a perturbative treatment of the 
disorder coupling 4 ^^ and no disordering transition as a 
function of the disorder strength is found. However, once 
topological defects are included in the analysis, the cou- 
pling of vortices to the random gauge field can lead to a 
disordered phase even at T = 0. This transition is driven 
by the creation of unpaired defects if the fluctuations of 
the gauge field are stronger than some critical value^ 4 ^ 
The critical disorder strength beyond which such defects 
appear can be estimated quite accurately when one cal- 
culates the free energy of an isolated defect in presence 
of disorder 4441 It turns out that a similar analysis of a 
single defect in a spiral in presence of disorder can be 
carried out with some modifications, at least at the level 
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of saddle point solutions. Within this approximation, the 
free energy of an isolated spiral defect is given by 

1 + (m z ) 2 b 

0F= ) \ lnR-[lnZ d ] D , (34) 

where the second term contains the corrections due to 
the disorder coupling, 

Z d = J ePycxp (-2 J d 2 x p k e ijk e abc n] n c k Q b ^j 

(35) 

with A^, n k obtained from Eqs. (|15ll . I|l(j|) an d 
(J20J. With use of the replica trick [lnZj^, = 
limjv^o jj In [^d] d> we nave ' assuming b < 0, 




with *„(x) = *(x - y„). We write 

J d 2 xd^ n d^ n , = -^A nn , + V 2 , (36) 

with V 2 ~ 27rlni? and A nn / ~ 47rln|y n — y n / For large 
separations |y„ — y n /| we approximate A„ n / ~ 47rlni? 
while for small distances A n „< is negligible. To find the 
highest weight configuration, the replicas are grouped to- 
gether in N / m sets containing each m replicas, with small 
distances between replicas within a set and large distance 
for replicas in different sets. [Z^] D then scales with R 
as 

[Z N ] ~ ^4A7rpj7V 2 +max m (2^--47rApiAf(Af-m)) 

In the limit N — > 0, maximization is replaced by min- 
imization with respect to m in the range < m < 1, 
so 

0F = [2 Pl ir - min < m <i (2/m + 4\p 2 irm)] InR . (38) 

For2\p 2 ir < lone finds (3F = 2[pi7r(l-2Api)-l] hi-R so 
that for pi7r(l — 2Xpi) < 1 free defects are favorable. This 
is the phase boundary for thermal creation of defects. At 
low temperatures, 2\p\-n > 1, one obtains j3F = 2irpi(l— 
a/8A/7t) InR and a critical disorder strength A c = 7r/8 
beyond which the disorder favors isolated defects even at 
T = 0. Similar considerations for the case b > lead 
to the same critical disorder strength and the condition 
7r(pi +P3)[1 — A(pi +P3)] < 2 for thermal creation of free 
defects. 

Let us first discuss the results for the disorder free case 
A = 0. The situation is summarized in Fig. [51 which 
shows the line separating the regime where free vortices 
exist from the regime in which all defects are bound. No- 
tice that the unbinding temperature goes linearly to zero 
in the limit b — > — 1. At b = — 1, free defects are present 



at any finite temperature. This is expected, as at b = —1 
and finite t s , the topological defects we discuss here lose 
their meaning as the stiffness for rotations around the 
collinear ordering axis disappears and the model becomes 
a 0(3)/ 0(2) model which has no finite temperature tran- 
sition. Whether or not free defects are present exactly at 
the point b = — 1, t s = depends on how this point is 
approached. To see this, we note that the symmetry of 
the model in the limit p^ — > 00 but finite p\ reduces to an 
XY symmetry as fluctuations of the 113 vector get sup- 
pressed which forces all fluctuations of the orthonormal 
pair ni 2 to lie within a plane. Therefore one obtains an 
XY model with stiffness p\ . In terms of the 6, t s parame- 
ters, this limit is approached as t s — > and 6^—1 with 
finite (1 + b)/t s — Ap\. Thus, depending on whether one 
approaches the point l + b = t s = with a slope larger or 
smaller than the critical one given by (1 + b)/t s = A/ir, 
one arrives at the disordered phase or the ordered phase 
of the XY model. This behavior is correctly reproduced 
by the free energy argument. The validity of the critical 
curve (1 + b)/t s = 4/ir also for finite 1 + b > is at least 
plausible, as topological defect solutions also survive in 
this limit. Below this line, the RG Eqs. I|20II21|1 hold and 
the system should scale towards the point 6 = 0. We can 
only speculate however what happens above that line. At 
least for some finite regime near b = — 1 the unbinding 
transition would presumably drive p\ to zero, as it does in 
the XY model, and affect the renormalization of p$ only 
weakly. Thus, the appearance of free defects will prob- 
ably modify the RG equations at high temperatures in 
such a way that the system will flow back to the collinear 
point b — — 1 as long as 1 + b remains small enough. For 
larger b the nature of the RG is unclear. Numerical sim- 
ulations on triangular Heisenberg models^S*^ have found 
however clear evidence for a defect unbinding transition. 
As the triangular Heisenberg model is believed to have 
initially b = 1^ it is likely that an unbinding transition 
indeed occurs for every initial value of b. As no RG equa- 
tions are available which can describe the transition, the 
form of the correlation length near this transition is un- 
known. It was however argued^ that the temperature 
dependence of the correlation length should cross over 
from the NLctM behavior to an XY behavior when the 
defects unbind. Numerical results seem to support such 

A 1 ? 

a scenario?' - 

Let us now turn to the case with disorder. Disorder 
will lead to the formation of free defects if A > tv/8. 
According to the free energy argument above, this critical 
disorder strength is independent of the stiffnesses p k and 
is thus also valid in the XY limit discussed above. For 
strong enough disorder, free topological defects will exist 
already at T = 0, invalidating our NLcrM analysis and 
producing very short low-temperature correlation lengths 
for the spiral. For XY models, the correlation length at 
T = behaves like £ oc exp(B/y/X — A c ) (with some 
constant B) near the critical disorder strength^ This 
form of the correlation length has a divergence of £ at 
A = A c which cannot be correct for the spiral because, as 
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FIG. 9: The critical line for the thermal unbinding of topo- 
logical defects is shown in b, t a space. 



discussed above, even without vortices, the coupling of 
any finite amount of disorder to the spins will lead to a 
finite correlation length. The correct dependence of the 
correlation length at T — on the disorder is expected 
to be an interpolation between the NLctM result and the 
XY behavior. 

Certainly, the free energy argument is not expected 
to work as well in the present model as it does for XY 
models. The parameters A and t s flow to strong coupling 
and thus the predictions of the free energy argument also 
become scale dependent. In other words, while at some 
small scale the system might look stable against the cre- 
ation of free defects, at some larger scale the system will 
become unstable according to the free energy argument. 
There does not seem to be a simple answer as to which 
scale is the correct one for applying the argument. Note 
that such problems do not arise in the XY model where 
the stiffness remains unchanged under the RG as long as 
vortices are ignored. In view of the divergence of the A 
and t s parameters in the NLtrM, one possible scenario 
would be that free defects will always be present at suffi- 
ciently large length scales. Numerical results do however 
not support such a scenario and rather point to the ex- 
istence of a finite critical temperature, 49 Below we shall 
apply the free energy argument with the bare parameters, 
i.e. at the smallest possible scale, which, if anything, 
would overestimate the stability of the system against 
free defect formation. 



IV. COMPARISON WITH EXPERIMENTS 

Let us now compare our results with experimental data 
on the SG phase of La2- x Sr x Cu04. Neutron scattering 
data— have revealed an incommensurability of the spins 
which scales roughly linearly with x. At very small x, a 
small deviation from the linear dependence is observed. 
Both features can be explained within the dipole model. 
The linear scaling is reproduced if the fraction of the 
dipoles which are ordered is doping independent, i.e., the 
number of ordered dipoles scales linearly with doping. 
The deviation from linearity might be explained with 
the increase of the average separation between dipoles 



at small x and a resulting diminished tendency of the 
dipoles to align. 

The same experimental data also shows the strong one 
dimensional character of the IC modulation, i.e. the in- 
commensurability is observed only in one diagonal of the 
Cu-lattice (6-direction) and thus breaks the symmetry 
of the square lattice. This phenomenon is usually inter- 
preted as being due to the existence of charge and spin 
stripes running along the other diagonal (a-direction). 
However, this IC is also expected for a spiral along the b- 
direction because its chirality breaks the translation sym- 
metry (it can spiral clock- or anticlock-wise) . In addition, 
this symmetry breaking is expected to show long-range 
order because the dipoles prefer a discrete set of lattice 
orientations. 

Another important consequence of the spiral chirality 
is the formation of topological defects. To judge, whether 
or not topological defects play a role in the LSCO SG 
phase, we need an estimate of A. We can use as a 
lower bound for A the result obtained from the collincar 
analysis- 4 where a disorder parameter equivalent to ours, 
but defined on the much smaller scale of the AF unit cell, 
was used. ^From a fit of the x dependence of the correla- 
tion length at x < 0.02 and large temperatures T > Tn, 
one obtains A ~ 20x. In this regime of x, the low tem- 
perature phase has long-range AF order and a collincar 
analysis is well justified. We assume that the linear de- 
pendence of the disorder parameter on x, A ~ 20x, also 
holds in the SG regime. This view is supported by mea- 
surements, which found that the width of the distribution 
of internal magnetic fields (i.e. local staggered moments) 
increases simply linearly with doping, with no detectable 
change on crossing the AF/SG phase boundary^ see also 
Fig. ^ It is remarkable that with our above estimate for 
the critical disorder strength A c = it/ 8 we find a crit- 
ical doping concentration x c ~ 0.02. Considering that 
A ~ 20x is a conservative lower bound of A at the long 
length scales relevant to spirals, we conclude that in the 
entire SG phase, free topological defects will be present 
already at T — 0, leading to a strongly disordered spiral 
phase. Experiments have in fact shown that the corre- 
lation lengths in the SG regime are rather short and of 
the same order as the periodicity of the IC modulation. 16 
While this is in accordance with the expected presence 
of topological defects, the correlation lengths are so short 
that the condition £ 3> |q s | _1 is not fulfilled. The regime 
where spiral correlations become dominant is therefore 
barely reached, and the RG scaling predictions cannot 
be well tested. 

While qualitatively the experimental data supports a 
description of the SG phase as a strongly disordered spi- 
ral state, both the extremely short correlation lengths 
and our limited understanding of topological defects pre- 
vent a more quantitative comparison. 

However, our suggestion that the incommensurability 
of the spins is related to ordered dipolar frustration cen- 
ters can be directly tested experimentally on co-doped 
samples La2_ x Sr x Zn z Cui_ z 04. Zn replaces Cu in the 
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CuC>2 planes and effectively removes one spin. Zn dop- 
ing leads therefore to a dilution of the AF but does not 
introduce frustration. Dilution is not very effective in de- 
stroying the AF order and pure Zn doping (with x = 0) 
leads to a destruction of long-range order only at perco- 
lation threshold that occurs for x « 41%^ Surprisingly, 
for very small Sr concentration x < 0.02 it was found that 
co-doping with Zn can increase T/v^i This is remarkable 
as both kinds of impurities lead to a reduction of Tjv in 
singly doped samples. A possible explanation for this 
behavior was suggested by Korenblit et alml They put 
forward an argument, that Zn impurities, if placed close 
enough to the localized hole state, will destroy the frus- 
trating nature of the hole bound state. While their micro- 
scopic picture of frustration is a classical one, a Zn impu- 
rity is also expected to strongly influence the properties of 
the bound hole state within a more realistic quantum me- 
chanical picture of frustration. Although Zn couples only 
weakly to the spin degrees of freedom, if placed near a Sr 
donor, it disturbs the symmetry around the Sr atom and 
modifies the nature of the bound hole state. As the Zn 
impurity breaks the sublattice pseudo-spin degeneracy of 
the bound hole, the orientation of the dipole moment is 
no longer annealed but becomes quenched. Another ef- 
fect of the breaking of the sublattice symmetry is that the 
weight of the bound hole wave function near the wavevec- 
tor (7r/2, 7r/2) or equivalent points will be reduced. As it 
is these wavevectors which are responsible for frustrating 
the spin background, one would expect a reduction or 
possibly a complete destruction of the frustration caused 
by the hole. Hence, the effective density of dipoles will be 
renormalized to x — > x(l— 72) where 7 must be calculated 
from a microscopic theory (experiments indicate that 7 
is of order 2)iS Co-doping with Zn then has two effects: 
First, it lowers the amount of frustration in the sample 
and thus increases the correlation length, which would 
explain the experimentally observed increase of Tjv with 
z for x — 0.017i£i*££ Furthermore, the effect of quench- 
ing the dipole moments will be the same as destroying 
them altogether with respect to the incommensurability, 
as the incommensurability is determined solely by the or- 
dered moments. Thus, co-doping with Zn will lead to a 
decrease of the incommensurability by a factor 1 — 72:. 
In contrast, within a stripe picture, co-doping with Zn 
is not expected to change the incommensurability as the 
hole density is not affected by Zn doping. Previous mea- 
surements in the superconducting phase [x = 0.12 and 
x = 0.14), where the stripe model is believed to be valid, 
have shown that the incommensurability indeed remains 
intact upon co-doping with Znm2£££^. Within a stripe 
picture, the only effect of Zn co-doping in the SG regime 
should be pinning of stripes, which would lead to a re- 
duced correlation length^ Therefore, neutron scattering 
experiments within the SG regime of Zn co-doped sam- 
ples could clarify the debate, if the magnetic incommen- 
surability observed in the SG regime is to be interpreted 
within a stripe or a frustration based model. 

It is interesting that symmetry arguments similar to 



those just used to discuss Zn co-doping also give a sim- 
ple explanation for the absence of any incommensurate 
signal in Li doped L^Cui-yLiyO^ For small y, these 
compounds show a magnetic phase diagram which is al- 
most identical to Sr doped samples^ with the notable 
exception that the magnetic correlations always remain 
commensurate^ Like Sr, each Li atom introduces an 
excess hole in the Cu02 plane which, at least for small 
doping concentrations, remains weakly localized to its 
dopant. The important difference is that Li replaces Cu 
in the crystal and thus has a different symmetry with re- 
spect to the magnetic sublattice ordering than a Sr hole. 
Specifically, the sublattice position of the Li atom breaks 
the pseudospin degeneracy present in Sr doped samples. 
Assuming that otherwise the origin of frustration is the 
same, the only difference between Sr and Li doped sam- 
ples is that the dipole moment assigned to the Li bound 
hole is quenched, whereas the one of the Sr hole is an- 
nealed. Thus, ordering of these moments and the devel- 
opment of incommensurate correlations cannot occur in 
Li doped samples. 

In conclusion, we have presented a detail picture of the 
dipole model of frustration and discussed its applicability 
to the weakly doped regime of cuprate materials. Most 
of the key characteristics of these materials were already 
known to be in accordance with the model and we showed 
that incommensurate correlations appear also naturally 
within the dipole picture. We extended the commen- 
surate model to allow for a description of the resulting 
disordered spiral spin phases. Finally, we suggested an 
experiment which would allow to verify whether the frus- 
tration based dipole model or the stripe picture is realized 
within the weakly doped regime of cuprates. 
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APPENDIX A: SU(2) REPRESENTATION 

The orthonormal basis n^. can be related to an element 
g of SU{2) through go k g~ x — n^- cr, or 

^^'^VV'} (Al) 
For the derivative one finds, using <9 M (gg^ 1 ) = 0, 

d»n% = ^tT{a a d ll ga k g- 1 +a a ga k d l ,g~ 1 } 

= \^{o k [g^o a g, g-%g]}. (A2) 
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Introducing g 1 d ll g = iA M - cr and with [<r J , cr- 7 ] 
2ieijk<J k one finds 



APPENDIX C: PROPAGATOR OF THE <p l 
FIELDS 



d„n% = 2e ijk A^. 



(A3) 



Therefore, we have (with p\^ = P2y) 



(A4) 



K + H A t 



with t 1 = 2(pi M + p 3/i ) and W 1 = 2(pi M - p 3/i ). 



APPENDIX B: EXPANDING THE ENERGY 
FUNCTIONAL IN tp* 

To do the RG, we introduce g — gexp(iip • cr ), where 
ip a are fast fields fluctuating with wavelengths [A _1 ,l] 
and g has only slow fluctuations in the range [0,A -1 ]. 
For the 1-loop calculation, we need to expand and 
up to second order in p a . We then find 

n" = ^tr {cr a gexp (i<p • a ) a 1 exp (— i(p • cr ) g^ 1 } 

= f% + ±i X {o*g[<p-tT,<x i ]g- 1 } 

+ jtr \v a g \y> ■ a a*<p ■ a -^(<p-crf a 1 



■0(p 3 ) 

= n? + 2e ijk ^h a k + ^^ h R% + 0{p% 



(Bl) 



where 



D«i ~ a 7. ( ~1 ~i ~M ^ „i „k i ^ „i „3 „k \ ~ — 1 



It turns out, that in the RG we will only need the diago- 
nal components of with j = k which have the much 

simpler form R™ — — 2 (e Z qi) 2 nf (we put here j = k = z 
to make clear that z is not a silent index, the equation 
also holds for j — x,y). Similarly, we find 



A' 



2i 
A 



tr {a 1 exp (-i<p ■ cr) [d^ + g l d ll g] exp (icp ■ cr)} 



„ , -tr j cr 4 (d^ip ■cr + — [ip-cr, d„<p ■ cr] 



+ i 



A^-cr,^ • cr 

2 



tp • cr A^-cr <p • cr 

1 r 



C% 3 ) 



= 1; + + e iik <pPd^ k + 2e l]k ^A\ - 2A\ <p 2 



As already mentioned, there is a small spatial 
anisotropy in the stiffnesses p kfJ ,, i-e. Pki + 1 Pki- We 
shall keep here the spatial dependence of the stiffnesses 
Pkfi up to first order in the anisotropy, assuming that 
the anisotropy k, which we define through p k i/Pk2 = 
1 + k, is independent of the k index. Thus we can ab- 
sorb the anisotropy into the i M parameter while b re- 
mains isotropic. We then define t s — \Jt\ti and i 12 — 
(1 ± n/2)t s . For future use, we also define the isotropic 
stiffnesses p k — -JPkiPk2- It is not clear whether the 
isotropy of b is preserved under the RG and we have made 
no attempt to write down the RG equations in presence 
of anisotropy. In principle, if b remains isotropic, the re- 
sults obtained below allow to determine the flow of the 
anisotropy parameter n under the RG. For possible fu- 
ture use, we will therefore keep the perturbative expan- 
sion with the anisotropy. The results used in the body 
of this work have however been obtained for an isotropic 
tfj,=t s , i.e. k = 0. 

We need to expand the exponential exp(—Hp) and in- 
tegrate out the p> % fields. Taking the average over the tp 1 
fields is done with the Gaussian term H v of Eq. (|18f) . The 
propagator for the p l is thus quite simple and becomes, 
to lowest order in the anisotropy k 



C*(x) := (^(xy(0)> c 



1 



k\ 



k 2 



2k 2 



d 2 k e ikx 
2(1 + bS iz )J (2tt) 2 k 2 

x(T(fc,A)-T(M)). (CI) 



The IR cutoff is provided by the function T(fc,A). A 
sharp cutoff, T(k, A) = 6(fc — A -1 ) has the disadvantage 
of producing a long-ranged C l and we therefore adopt 
instead Y(fc,A) = [1 + (fcA) -2 ] -1 , which renders C l short 
ranged. 

In our RG calculation we will mainly need C 4 (0) which 
has the form 

C x (0) = C*(0) = — In A+C(k 2 ), C z (0) = — -6^(0). 

47T 1+0 

Another useful formula is 



t~ x j d 2 x (d,C*) 2 = l -C*{0)+O(K 2 ). 



APPENDIX D: RENORMALIZATION 



(C2) 



+ 2A M - ¥ > i P i +0(<p 3 



(B2) 



We can immediately discard all terms of third or higher 
power in A^ as these terms are irrelevant in a RG sense. 
Terms second order in A M renormalize t„ and 6, whereas 
terms linear in A M are responsible for the renormalization 
of the disorder variance A. 

First, we note that the terms H 2 , H 3 do not contribute 
to the renormalization, as was pointed out for the calcu- 
lation of the RG for the disorder free system in Ref. [ 43 ]. 
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This is because these terms are linear in cp while they do 
not involve a disorder field Q M . For an abelian theory, 
such terms cannot contribute because the fast <p % fields 
and the slow A M fields have their support in orthogonal 
parts of the wave vector space. Here, for the non-abelian 
case, this argument is not sufficient because the A M fields 
are not linearly related to the fields g. For the present 
non-abelian theory this is nonetheless true, although an 
explicit calculation is required to see this. For example, 
iff does not contribute, because its contribution is built 
from terms of the form (we omit the upper i indices of 
C % and A 1 ^ here for simplicity) 

J d 2 x J dW i M (x)i^(x')a M ^C(x - x') (Dl) 

To evaluate this term, we change to center of mass (y) 
and relative (y') coordinates and then perform a gradient 
expansion in the relative coordinate. Only the lowest 
order contribution is of interest, as higher order terms 
involve a local coupling of the type (d v ) n A^i with 
n > which are irrelevant from a scaling point of view. 



The lowest order term is then 

- J d 2 y 4(y)4,(y) J d 2 y'd^,C(y') (D2) 

which vanishes because the last integral is zero. In the 
following we will omit H 2 and H 3 from the analysis, be- 
cause terms involving them do not contribute. This can 
be shown for each term in a way similar to the one just 
shown. 

We want to find the RG equations up to second order 
in tu and A. In the nth order of the cumulant expansion 
of j r , Eq. (|19fl . we only need to consider terms which 
have a total number of cp and Q p fields less than 2n + 2. 
This is because each term of order n carries a factors t~ n 
from the prefactors of the terms in H p and each pair of 
<P (Q/j) produces a factor t s (A). 

We begin first with the terms renormalizing and 6, 
where we give a detailed calculation only for the terms 
up to second order in H p . The calculation of higher order 
terms is quite lengthy although conceptually easy and we 
therefore just present the results of the calculation. 



I 

1. Terms which renormalize t M and b 



a. First order in H p 



There is only one term quadratic in A p which contributes, H4 (the tp l average over H3 is zero). 



(H 4 ),„„ = -4— / d 2 x 



= -4- 



{e zjk f(A k A C"'(0)-(I* £C(0)+ (A*) C*(0) 



1* 1* 
2 



-46 1: 



d 2 x 



C x (0). 



(D3) 



b. Second order in H v 



Terms with odd numbers of cp 1 or Q M are zero after performing the ip 1 and disorder average. There are then only 
two terms we need to consider, H\ and i? 2 x (ff 2 3 has a total of six ip 1 and Q 1 ^ fields and does not contribute and H 2 
terms do not contribute as mentioned above). For H\ we have 



1 r 



(Hi) 



tpc 



D 



L (Hi) 



z > V 



(D4) 



<i 2 x d 2 x' AJ,(x)i4*,(x')ey*ei/j/v (1 - bS lz + 2b6 jz ) 



x (1 - U Vz + 2bSj' Z ) (^^'(x)^(x)^^'(xO/'(xO) 



The four point average can be decomposed according to Wick's Theorem. Nonzero contributions arise from the 
contractions (jk r ) (j'k) and (jf) (kk'}. We again employ an expansion of i? 2 in the relative coordinate and keep only 
the zeroth order term of the expansion. This yields 



2 ( H D V c 



2t: 



d 2 x.A l ll A l tl e ljk ei>j>k> (1 



+ 2bS jz ) (1 



2bS fz ) 
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x (SjpSkk, - S kjl 5 jkl ) / d 2 y d^C\y)d^C k (y) 



4t~ 



d 2 x 4 (e ijk f(l 



+ 2bS jz ) (1 - b5 lz + bS ]Z + b8 kz ) 



x / d 2 y d^(y)d ll C k (y). 



With use of Eq. ljC2)l . we finally find 

The other second order contribution is 
1 



2t~ x / d 2 x 



A*Q. + b) + [Al) 6(6-3) 



C x (0) 



8 J d 2 xd 2 x' Pk^Pk' fj,' tijkti' j' k< ^abc^a'b' d { 



£kimn°;n c m A 



+e nrm ,ni,ni,A™:)5 w C ; (x-x') [q* (x)Q»',(x') 



Using 



Q£(x)Q£,(x') 



(D5) 



(D6) 



(D7) 



find after some algebra 



<5&f/ <W <5(x- x')A, eabcta'bd = Saa> 5 CC > - 5 ac >8 ca > and the orthonormality of the n fc , we 
1 



(Hi,) 



= 2A Y t: 



D 



A: 



C x (0). 



(D8) 



Higher order terms can be evaluated in much the same way as the first and second order terms, although the large 
number of indices makes their evaluation more tedious. We therefore refrain here from a detailed presentation of these 
terms and just state the results. 



Third order in H v 



Terms of second order in A 2 are produced by (Hi + H cl + H c3 ) (H c2 + H c4 ). However, only the terms Hi(H c i + 
H c z)(H c2 +H c4 ) have even powers of Q M . Terms with eight or more <p and Q M fields again do not contribute to second 
order in A,i M . Thus we are left with only H\H c \H c2 . We find 



\H\H c iH C 2 



-2\t~ 2 b 

h 1 



d 2 x 



',(! + &) + (l*) 2 (6 -3) 



C x (0). 



(D9) 



We further need to consider terms of the type (H c2 + H c4 ) 2 H 4 . Only H 2 2 H 4 has less than eight </?,Q/j fields and even 
powers of both fields. We find 



1 



(H 2 2 H 4 ) 



2\ b t- 1 t- 1 



d 2 x 



3 [A 



C x (0). 



(D10) 



Fourth order in H v 



Possible contributions arise from the terms (Hi + H c \ + H c ^) 2 (H c2 + H c4 ) 2 . Discarding terms with ten or more 
(p l ,Qfj, fields, we are left with H 2 2 H\ and H 2 2 H 2 X . However, the connected part of the ip l average of H 2 2 H 2 X is zero (its 
finite disconnected parts enter the renormalization of the disorder, see below), and the only contribution is therefore 



H c2 H 



= AC 1 *" 1 



d 2 x 



A*(2 + 6)(l + 6) + (A*) 6(6-7) 



C x (0). 



(Dll) 



Terms of the form H 4 (H c2 + H c4 ) 3 do not contribute because their disorder average is zero. Higher order terms in H p 
do not contribute because they either involve more than four Q M terms and are therefore of higher order than A 2 or 
they do not contain finite connected parts. For example, the term (H 4 H^ 2 ) r decomposes into products of averages 

of (H 4 ) vc or (H 4 H 2 c2 ) vc and (H 2 c2 ) ^ ^ 



' tpc 



Terms which renormalize A 



linear in A 1 . We list the contributions order by order 



To find the renormalization of the variance of the dis- 
order distribution, we first collect all connected terms 
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below. 



c. Third order in H„ 



a. First order in H p 

Only three terms are linear in A 1 , H\, H c \ and H c3 . 
However, both Hi and H c i have a zero ip z average and 
only (H c3 ) vc contributes. 



There are contributions from {H c iH c2 H c4 ) , 
(H c3 H? 2 ) vc and (H^)^. The terms (H c3 H? 4 ) vc 
and (^H\H^ c do not contribute, as they contain eight 
or more Q^, <p l fields. 



b. Second order in H„ 



d. Fourth order in H„ 



At second order there are contributions from 0nl y onc term contributes, (H x Hl 2 H cA ) ^. All other 
H c iH C 2) and (HiH C 4,) . There is no contribution to terms have ten or more Q^, ip* fields or more than three 
second order in A, t M of the disorder renormalization from Q M fields and thus do not contribute. The same argument 



(H c3 H ci ) vc because this term has six Q 1 ^, tp l . 



applies to all terms generated by higher order of H p . 



J 



3. Calculating the renormalized disorder variance 



We now must calculate the variance of all terms at the new length scale A 1 which are linear in A^. These are the 
terms just found above plus H c0 . Thus, we need to calculate the variance of 



-H c0 - (H c3 ) vc + (H cl H c2 ) vc + {H x H c ^ vc - (H cl H c2 H c4 ) vc - - (H^ 2 H c3 )^ 



~2 ( H c2 H l) vc + 2 ( H l H c2 H c4) vc 

To order A 2 , the following terms contribute to the variance. 



[ H co] D 



\t~ z I <f x 



A] 



(i*) 2 (l + 6) 2 }, 



D 



d 2 x 



(o 2 +(^y 

- (Al) 2 26(1 + 6) lc*(0), 



(H 1 H c4 ) vc H c0 ] d = -4A^ 2 |d 2 x|[(i-) 2 +(^) 



(1 + 6) 



+ 



(i*) 2 (l-6) 2 (l + 6)W(0), 



(H c1 H c2 H c4 ) vc H c0 ] d = 2AV|d 2 x|[(^) 2 +(i^) 



6(1 + 6) 



A') 26(6 2 - 1) C*(0) 



(Hl 2 H c3 )^ c H c0 



D 



(^) 2 26(l + 6)|^(0), 



(H 1 H 2 c2 H c4 ) vc H c0 ] d = -2AVt s - 1 |d 2 x|[(i-) 2 +(l^ 



(l + 6)(2 + 6) 



(A;) 2 2(1 + 6)(l-6) 2 }c*(0) 



{H cl H c2 } 



d 2 x 



A' 



(D12) 
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+ {A^ 2 (2 + t s t- 1 )b 2 \c x (o), 



( H ^)l] D = A a t-V/^{[(^) a + (A») 



+ 



(i*) 2 (l-6) 2 W(0) 



(H 1 H^) ipc (H cl H ta ) vc ] D = -2AV^ 1 /^ 2 x{[(^) 2 + (i)l) 



6(1 + 6) 



+ (i^) 2 26(6-l)|^(0). 



The sum of the above terms is (we now again set = t s 

2 / - \ 2' 



1 



At7 2 d 2 x 



a: 



1 + M uMa 



+ (i-) 2 ((l + 6) 2 - 



2 4(l + 6) 3 t s + (3 + 6 6 + 6 2 )A cx(0) 



(D13) 



a. On t/ie calculation of disorder terms 
As an illustration, we give details for the calculation of the variance terms for a relatively simple term, 



and a more involved one, 



(h c1 h 2 2 h c4 ) vc h, 



rO 



. For 



D 



c() 



we have 



D 



d = 8 j d 2 ^d 2 ^p k ^p k >^e tjk ee^k'^abc^a'Vc'C\G)hfni,A^ x 

\ c\ ~ a ~ c ~Xi , c\ ~ c ~ a 7 m , r. ~ a ~ c 2m 

< z€ji 7n tkiqn m n q si^ -f- ztii 7n €ki q n q nj -+- z,eii m tjiqn q n k /±^ 
-n^ngAj, ((e i(m ) 2 + (e,/ m ) 2 + (6 Hm ) 2 )} [Q»(x)Q^(x')] D 

{ (^jimekiqh^h^A^ + e i;m e fei9 ^n^i™ + e^e^n^i") C'(0) 

-fi«n^4 (2 + (l + 6)- 1 )^(0)} 

= 16A y d 2 x(e yfc ) 2 (4) 2 { W ^C 4 (x)+p 2 At C 4 (x)+ Wv ^(x) 

+Pfc M 7VC fe (x) +p^C fe (x) +p fe ^ AI C fe (x) 
-(2 + (l + 6)- 1 )C a; (0)( ? , 2 M + ? , fe ^)}, 



(D14) 



where we again used the orthonormality of the rife. Performing the summation over the silent indices, one finally 
obtains 



= 16A J d 2 x| (i* 
= 4At; 2 | d 2 



A l 



{Pi Pi) + (K) 2 ( 4 Wim - 4pU } ^(0) 
i^) 2 + (i^) 2 ] b (i*) 2 26(1 + 6) J C*(0). 

We now turn to the more lengthy evaluation of (_ffi_ff 2 2 _ff c4 )^ c H c0 . We have 

HiH 2 2 H C 4 = 16 y <i 2 X(i 2 x'd 2 x''(i 2 x'''pfe M t^/pfe/p/pfe//p//pfe/// M ///eijfeei'j/fe'ei''j//fe//ei///j///fe//^ x 
eo6ceo'6'c'eo"6''c"eo'"6'"c"'^' (1 - + 2bSj< z ) h*,, h c k ,,hj,,, h c k ,,, x 
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(D15) 



We now need to perform the average over the ip fields. For convenience, we split HiH^Hd — A + B into two terms, 
where A corresponds to the part of HiH^ 2 H c i which involves the first term in the curly brackets in Eq. I|D15|) and B 
corresponds to the second term in the curly brackets. For (A) , we need to calculate the average 



(s M ^(x)/(x)a M v''(x')^ fc '(x')V^''(^')V'/''( x/// )) 



(D16) 



which can be easily done via Wick's Theorem. However, not all possible permutations of pairings will contribute. 
All terms involving either of the contractions (id) or (j'k 1 ) vanish as d^C x (0) = 0. Although not immediately 
apparent, terms involving the pairing also do not contribute to one loop order. This can be seen only after 



the computation of the disorder average 



and a gradient expansion similar to the one employed below 



Eq. (|D1|I . Using the same arguments as we used for the term I|D1|) . all contractions can then be shown to give 

no contribution. Furthermore, all contractions which are identical up to a permutation of the indices i" and i'" will 
give the same contributions after the disorder average is taken, as discussed below. We therefore only write down half 
of the permutations and indicate the others by {// <-> ///}. Thus, we only need to keep the following terms, 

(9 M ^(x)^(x)^^'(x')¥' fe '(x')^»/'(x")^»/"(x"'))^ 

6 ii ,8w8 k , v „d IM d llf C i {x. - xX-C d (x _ x ")^»<C < "'(x' - O 
+6 ik ,5 di „5 fi ,„dt l C i (x - x')9 M "C d (x - x")c^^C J "'(x' - x '") 
+6 u »6 d:l ,5 kH „,d^„C i (x - x")^C d (x - xOVC"'^ ~ x"') 
+5 u „5 dk ,8 j , i ,„d^,C i {x - x")C d (x - x')0 /1 ,0 /t „,C 4 "V - x"') 



+ {// «-> ///} 

Let us now perform the disorder average (A) vc H c o 



Q*(x)Q»(x)Q£,,(x")Q£„,(x'") 



For this, we need to calculate 



(D17) 



(D18) 



where the variables carrying a tilde arise from the H c q term. Again, we can use Wick's Theorem to decompose the 
average. Of the three possible permutations of pairings, two involve either of the two contractions (66") or (bb 1 "). 
Neither permutation contributes. This is easily seen for the (bb") contraction and the explicitly written terms in 
(|D17() because they all involve after the contraction a derivative of C x (0) and thus vanish. The same terms also 
do not contribute for the case of a (bb'") contraction, which again can be seen with a gradient expansion and using 
arguments analogous to those below Eq. I|D1|) . Therefore, only one term of the disorder average must be kept, 



Q|(x)Q*(x)Q*:(x'0QOx''') 



D 



A 2 <5^ W<W<W"<Kx - x)<5(x" - x'") 



(D19) 



The terms in (|D17(1 which only differ by a permutation of the double primed and triple primed variables give then 
identical contributions, as such a permutation simply relabels the variables associated with the two H C 2 terms in 

(A) H cQ ] . With llDT7llDT9)l we then have 



D 



128A 2 

(1 — bv t ' z T ^wuj'z 

x 

^(W'W^C^x - xX"C d (x _ x")^»C 4 "(x' - x") 

+6 ik ,8 di „5 fi „d l ,C i (x ~ x')^»C d (x - x")a M / VC*"(x / - x") 
+6 u »6 d:l ,6 k , l „d^,C l (x - x")^C d (x - x')d^C l "(x' - x") 
+6 ii »6 dk ,5 j , i „d fl d fl „C i (x - x")C d (x - ^(/(j' - x") 



(D20) 
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The integration over x" can now be performed with 

t~}, J #x!%r,C°(x - x'X"C x (x' - x") = ic a (x - x') + 0(k' 



(D21) 



The remaining double integral over x and x' can then again be approximated with a gradient expansion in the relative 
coordinate and employing Eq. (|D21|) . We then obtain (we denote the center of mass coordinate again by x) 



16X 2 t s I d 2 x (p\,, + p k "Pj") (ei"j"fc») 2 ti'j'k'A^ (1 - b5i, z + 2b5j, x ) x 

0i/3 d e idk {AfypiuPkn + A*p dft p kfi + A^p 2 ^ + A^p^pa^ x 

[Sij'5di"Sk'i" Pk' — Sik'8di"8j>i"(3j' — 8u"5dj'Sk'i"(3k' + §u"fidk'fij'i"Pj'] x 
xC x (0). 



(D22) 



where Pk is defined through Pi = Pi = 1, Pz = (1 + b) 1 and Pftt s /t^ — Pku,- After some straightforward algebra, one 
finally finds 



(A)H c0 



D 



32X 2 t~ 2 t 3 s I d 2 x 



((pi + P3) 2 + ((pi + P 3 ) 2 + + 2 w*) + 

(A*) 2 8(1 - 6) (p? + p lP3 ) (p! + P3 A C*{0). 



The calculation of 



cO 



can be done in much the same way as just shown for 



D 



-32X 2 t~ 2 tl / d 2 



(D23) 
One arrives at 



D 



fl s 



+ (Al)\{l-b)p 2 (p 1+ pzf\c*{G) 



(Pi + P3) 2 + "TTT ) (Pi + P3) 2 



1 + 6 



(D24) 



Finally expressing all p k through b and i s , one obtains for 
(H 1 H 2 2 H c4 ) H c0 



(A + B) H c < 



D 



D 



(^) + ( A l) 
+ [A^Y 2(1 + &)(!- bAc(0). 



(l + 6)(2 + 6) 



(D25) 
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